Charge-density-waves and superconductivity as an alternative to phase separation in 

the infinite-U Hubbard-Holstein model 



F. Becca, M. Tarquini, M. Grilli, and C. Di Castro 
Dipartimento di Fisica and Istituto Nazionale per la Fisica della Materia 
Universitd La Sapienza, 
Piazzale Aldo Moro, 00185 Roma, Italy 



o\ : 
u ; 

Oh 
<\ 

in ■ 



> 
in 
m 
O 

o 

ON 
-I— > 



o 
o 



X 



We investigate the density instabilities present in the 
infinite-U Hubbard-Holstein model both at zero and finite mo- 
menta as well as the occurrence of Cooper instabilities with a 
specific emphasis on the role of long-range Coulomb forces. In 
carrying out this analysis a special attention is devoted to the 
effects of the strong local e-e interaction on the e-ph coupling 
and particularly to both the static and dynamic screening 
processes dressing this coupling. We also clarify under which 
conditions in strongly correlated electron systems a weak ad- 
ditional interaction, e.g. a phonon-mediated attraction, can 
give rise to a charge instability. In the presence of long-range 
Coulomb forces, the frustrated phase separation leads to the 
formation of incommensurate charge density waves. These in- 
stabilities, in turn, lead to strong residual scattering processes 
between quasiparticles and to superconductivity, thus provid- 
ing an interesting clue to the interpretation of the physics of 
the copper oxides. 



I. INTRODUCTION 

Besides the large critical temperatures, the supercon- 
ducting copper [Oxides display many anomalous normal- 
state properties!!!. The understanding of these properties 
is not only a fascinating theoretical challenge, but would 
also shed light on the pairing mechanism leading to high 
temperature superconductivity. 

The anomalous properties of the normal phase have 
been interpreted along two distinct theoretical lines. The 
low dimensionality of these highly anisotropic systems 
and their correlated nature have been proposed to be at 
the origin of a breakdown of the Fermi liquid (FL) . In paij-i 
ticular the concept of Luttinger liquid in two dimcnsioncl 
was put forward as a new paradigm for the normal state 
of the copper oxides and it was intensively investigated!! 
However, it was recently showrO that the Luttinger liq- 
uid is only stable in one dimension. Above one di- 
mension, the Fermi liquid picture is recovered when the 
bare electron-electron (e-e) interaction is non singular. 
This result would support the alternative aptitude, which 
has been to accept the Landau theory of normal FL's 
as a suitable starting point. The anomalous proper- 
ties would then arise as a consequence of singular scat- 
tering processes at low energy between the quasiparti- 
cles. Along this line magnetic scattering has been con- 
sidered to be responsible for both the anomalous prop- 



erties of the normal phase and for the superconducting 
pairing!]. It was also proposed that excitonic. scattering 
could give rise to the so called marginal FLH, and pro- 
vide a pairing mechanism. Singular scattering is also 
obtained by gauge fieldsu, which arise by implementing 
the resonating-valence-bond idea in the t-J model. 

The above theoretical lines have more recently been 
joined by a different scenario suggesting phase separa- 
tion (PS) as a possible source of anomalous scattering 
and, therefore, of anomalous normal-state behaviorBc. 
Emery and KivelsonB suggested that, although long- 
range Coulomb (LRC) forces spoil PS as a static ther- 
modynamic phenomenon, the frustrated tendency to- 
wards PS may still be important and give rise to large- 
amplitude collective density fluctuations. Approaching 
the problem within a coarse-grained model, they sug- 
gested that these fluctuations may be responsible for the 
anomalous behaviour of the normal phaae and for the 
superconducting pairing. In a recent workH, two of us as- 
sessed the relevance of charge instabilities (PS or charge 
density waves) as a mechanism for anomalous scatter- 
ing, by determining the dynamical effective scattering 
interactions among Fermi-liquid quasiparticles close to a 
charge instability, both in the presence and in the ab- 
sence of LRC forces. This analysis consisted in a mi- 
croscopic treatment of the Hubbard-Holstein model in 
the infinite-U limit, finding that, both in the presence 
and in the absence of LRC forces, the dynamic effec- 
tive interaction have a singular behaviour, strongly af- 
fecting the single-particle and the transport scattering 
time. This scenario is obviously sensible if i) the consid- 
ered microscopic model displays PS for some parameter 
region and ii) the real copper-oxide systems actually are 
in the proximity of a charge instability!!!]. As far as point 
i) is concerned, PS seems to be a rather generic and ro- 
bust phenomenon in the context of strongly interacting 
systemstll. Indeed, after PS was shewn-to be present in 
the phase diagram of the t-J mode£3il3, it was pointed 
out that PS-C Qrampr ily occurs in models with short range 
interactionll3~t30l23 provided the strong local e-e repul- 
sion inhibits the stabilizing role of the kinetic energy. 
Moreover it was repeatedly claimed that PS and super- 
conductivity can be related phenomena-irrespective of 
the nature of the short-range interaction!!!!. 

On the other hand, the frequent occurrence of PS in 
models of strongly interacting electrons is made intrigu- 
ing by the observation of PS in oxygen-doped supercop. 
ducting copper oxides (La2Cu04+ y ) of the 214 classl!2l 
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Although the electronic origin of PS in these cuprates 
is still to be established, the contemporary presence of 
a strong e-e interaction and of PS in a real system and 
the robustness of the PS concept in theoretical models 
is suggestive. The reason why only La2Cu04+ y seems 
to phase separate is that LRC forces effectively oppose 
the separation of charged particles. Only when the nega- 
tively charged oxygen ions are sufficiently mobile the pos- 
itive holes can separate being accompanied by the oxygen 
countercharges which compensate for the charge unbal- 
ance. Nevertheless, even in those systems where LRC 
forces are present to prevent a thermodynamic instabil- 
ity, phase separation may remain in the system in the 
form of a tendency toward charge, ags^g^t ion, possibly 
giving rise to supercoiad^ctivitju'E^rtBtilcJ or to anoma- 
lous normal propertiestffl. In particular it might well hap- 
pen that the long-wavelength density fluctuations associ- 
ated to PS are suppressed in favor of shorter-wavelength 
density fluctuations giving rise either to dynamical slow 
density modeaa. or to incommensurate charge density 
waves (CDW)til. This latter possibility was recently 
put forward to explain neutr«Hi-scattering results in a 
Lai.4sNdo.4Sr .i2Cu04 sampleE2l. In this case it was 
proposed that the (low-temperature tetragonal) lattice 
structure and the filling (close to 1+1/8 holes per Cu02 
cell) were suited to pin the density fluctuations giving rise 
to a static CDW phase. The formation of striped pat- 
terns in the CuC>2 planes of Bi221|2,compounds was also 
shown from EXAFS experimentscJ Local density fluc- 
tuations could also account for some r esult s of neutron 
scattering experiments in 123 materialsESc3 

Based on the observation of an antiferromagnetic phase 
close to the superconducting one in the phase diagram of 
the copper oxides, previous analyses put emphasis on the 
role of magnetic-, coupling in originating the slow den- 
sity fluctuations^. However the generic occurrence of 
PS in theoretical models with different interactions, in- 
dicates that a definite choice of the mechanism leading 
to PS could be misleading or premature in the absence 
of more stringent experimental indications. Moreover, 
while in the models mentioned above the additional in- 
teractions inducing—PS are of purely electronic origin, it 
was shown in Ref.Ell that, in the presence of a strong 
local repulsion, also the lattice may introduce an effec- 
tive attraction determining PS in the Fermi liquid. This 
latter result showed within an infinite-^/ three-band Hub- 
bard model that the instability occurred for reasonable 
values of the e-ph coupling indicating that PS by no 
means requires unlikely parameters, unusual mechanisms 
or purely electronic interactions, but it can simply re- 
sult from the interplay between the strong local repul- 
sion and the (weak) additional attraction provided by 
the lattice. This theoretical observation is accompanied 
by some experimental evidences that the lattice can play 
a non-negligible role in determining the physics of the 
superconducting cupratesEH. In particular a sizable cou- 
pling between some phonons and the carriers is implied 
by the presence of polaronic effectsES for the very lightly 



doped compounds, by the copper and oxygen isotopic ef- 
fect present in La2- x Sr x Cu04, by the Fano line shapes 
in Raman spectra and by the rather large frequency shift 
and linewidth broadening of some phonons at T c . 

In the present paper we-tmrsue the investigation along 
the route opened in Ref.cJ The occurrence of a phase- 
separation instability was justifyied within a general 
Fermi liquid analysis demostrating that the strong in- 
teraction is responsible for vertex corrections, which are 
strongly dependent on the VFq/u) ratio, where vf is 
the Fermi velocity and q and u> are the transferred mo- 
mentum and frequency respectively. These corrections 
generically lead to a strong suppression of the effective 
coupling between quasiparticles mediated by (a single) 
phonon exchange in the VFq/u) 1 limit. However, such 
effect is not present when VFq/w <C 1, which is the rele- 
vant limit for the effective interactions entering the Fermi 
liquid expression for the compressibility. In these effec- 
tive interactions in the dynamical limit, the e-ph coupling 
is therefore not effectively screened opening the way to a 
possible violation of the stability criterion for the Fermi 
liquid in some regions of the parameter space. In Ref.cJ 
a detailed analysis was then carried out using a slave- 
boson approach for the infinite-U three-band Hubbard 
model describing the basic structure of a Cu02 plane in 
copper oxides. In the presence of a coupling between the 
local hole density and a dispersionless optical phonon, 
it was explicitly confirmed the strong dependence of the 
hole-phonon coupling on the transferred momentum ver- 
sus frequency ratio and it was also found that the ex- 
change of phonons leads to an unstable phase with nega- 
tive compressibility already at rather small values of the 
bare hole-phonon coupling. Close to the unstable region, 
Cooper instabilities both in s- and d-wave channels were 
detected supporting a possible connection between phase 
separation and superconductivity in strongly correlated 
systems. 

We now start from the infinite-U single-band Hubbard 
model in the presence of,-an optical phonon coupled to 
the local electron densitycS. Due to its relative simplic- 
ity with respect to the three-band Hubbard model we will 
be able to extend the model in a rather direct and straigth- 
forward way so as to include the LRC forces between the 
electrons. This extension is particularly important since, 
as mentioned above, LRC forces obviously affect the oc- 
currence of PS instabilities and could provide a clue in 
explaining the relative rarity of this phenomenon in the 
real materials. In this way, as briefly reported in Ref.Q, 
we also provide a microscopic derivation of an incommen- 
surate CDW instability directly from a system of strongly 
correlated electrons with all the physical implications in- 
dicated above. Therefore this topic represents a key issue 
of our investigation and may definitely be considered as 
the main point of our analysis. 

In Section II we introduce the model and the formal- 
ism. Readers, who are not interested in technical details 
can directly move to Section III, where we present the 
results concerning the physical properties of the model 
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in the absence of LRC forces. The effects of LRC inter- 
actions are reported in Section IV, which thus represents 
the core of the present paper, while in Section V we dis- 
cuss the results and we draw our conclusions. 



II. THE HUBBARD-HOLSTEIN MODEL 

Our starting point Hamiltonian is the twodimensional 
Hubbard model with an additional dispersionless phonon 
mode A coupled a la Holstein 

H = ~t ]T (c] CT c iCT + H.c) -t> ]T (clc ja + H.c 



At the mean-field [N = oo) level, the model of Eqs.(g)- 
(0) is equivalent to the standard, purely electronic U — 
oo Hubbard model without coupling to the phonons, 
which has been widely considered in the literatures. 
In fact, at mean field level no role is played by the 
phonons because our electron-lattice coupling depends on 
the difference between the local and the average density 
and this difference naturally vanishes in the mean-field 
approximation^. The average number of particle per 
cell is n = (1 — S)N/2 and 5 = corresponds to half- 
filling, when one half electron per cell and per spin flavor 
is present in the system. 

The mean-field self-consistency equations are obtained 
by requiring the stationarity of the mean-field free energy 
and they determine the values of b\ = Nrg = (bi) 2 and 
of Aq = (A;). Then the mean-field Hamiltonian reads 



u 



£ A]A t + g ]T (4 + Ai) (n% a ~ <rw» , (1) 



H M f = Eu(y ka Cka - (mo - A ) c\ a cka 



where and indicate nearest-neighbor and 

next-nearest-neighbor sites respectively and — c\ a Ci G 
is the local electron density. Since we are interested in the 
limit of strong local repulsion we take the limit U — > oo, 
which gives rise to the local constraint of no double oc- 
cupation rii a < 1. To implement-this constraint we 
use a standard slave-boson techniqueE3-Ej by performing 
the usual substitution c\ a — j—pj^fa;, a a —* b\a a . We 
also use a large-iV expansion^ in order to introduce a 
small parameter allowing for a systematic perturbative 
expansion without any assumption on the smallness of 
any physical quantity. Within the large- N scheme, the 
spin index runs from 1 to N and the constraint assumes 
the form J2 a c L c i<r+&I&i = y ■ A suitable rescaling of the 
hoppings t — > t/N and t' — ► t' /N must, in this model, 
be joined by the similar rescaling of the e-ph coupling 
g — > g/VN in order to compensate for the presence of 
N fermionic degrees of freedom. The model can then be 
represented as a functional integral 
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- » £ + A\) (n la - (n la )) + u; £ A^. (4) 

i.a i 

where a local Lagrange multiplier field A^ has been intro- 
duced to implement the local constraint. 



(5) 



where E k = —2tr\ek is the quasiparticle band with 
£fe = (cos k x + cos k y ) + a (cos(fc x + k y ) + cos(k x — k y )) 
(we define a = t'/t). In particular it turns out that the 
square of the mean-field value of the slave-boson field 
bo, 6q = Ntq = NS/2, multiplicatively reduces the hop- 
ping, t — > tb 2 ,, thus enhancing the effective mass of the 
quasiparticles. Moreover, at this level the single parti- 
cle self-energy does not introduce a finite quasiparticle 
lifetime. Then, in this model the single-particle Green 
function of the physical fermions at N = 2 has a quasi- 
particle pole with a finite residue given by the square of 
the mean-field value of the slave-boson field b\ . Thus for 
any finite doping S the system is a T=0 Fermi-liquid. 

On the other hand, at half-filling bo = S = and the 
system is insulating with a vanishing value of both the 
quasiparticle bandwidth (i.e. an infinite quasiparticle ef- 
fective mass m*) and a vanishing residuum of the polar 
part in the single-particle Green function. 

As far as Ao is concerned, this quantity rigidly shifts 
in a doping dependent way the bare chemical potential 
Ho and is self-consistently determined by the following 
equation 

\o = \ o + a\ 1 = 2tJ2.f(Ek)£k 

k 

= 2tY J f{E k ){h + a lk ) (6) 
k 

where f(E) is the Fermi function and (3k = cos k x +COS k y 
and = cos(k x + k y ) + cos(k x — k y ). 

The presence of the coupling with the phonons intro- 
duces new physical effects when one considers the fluc- 
tuations of the bosonic fields. Since only a particular 
combination a = (A^ +A)/(2yN) of the phonon fields A 
and A^ is coupled to the fermions, it is more natural to 
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use the field a and to integrate out the orthogonal combi- 
nation a = (A — A')/(2VN). Then the quadratic action 
for the boson field a reads 

i? Ph o„ = ^E^ ± ^ a l a - ( 7 ) 

— ^0 

where we have transformed the imaginary time into 
Matsubara frequencies. .-Moreover, it is convenient to 
work in the radial gaugeO, the phase of the field hi = 
y/Nri exp(— is gauged away and only the modulus 
field Ti is kept, while Xi acquires a time dependence 
\ — ► A.j + d T 4>i. Thus one can define a three-component 
field A^ = (Sr, 5X, a) where the time- and space- 
dependent components are the fluctuating part of the 
boson fields ~ tq (1 + 5r{), Xi = —iXo + SXi and dj. 

Writing the Hamiltonian of coupled fcrmions and 
bosons as H = Hmf + H\, os + H- lat , where Hmf is the 
above mean-field Hamiltonian, which is quadratic in the 
fermionic fields, -ffbos is the purely bosonic part, also in- 
cluding the terms with the a, r and A bosons appearing 
in the action (||) and in H p \ lou , Eq.(Q). iJi nt contains 
the fermion-boson interaction terms. The single-band 
U = oo Hubbard model also contains a four-leg vertex 
arising from the hopping part of the Hamiltonian (see the 
second term on the r.h.s. of Eq.([|)). The two fermionic 
legs of this vertex can be contracted (see Fig. 1) giv- 
ing rise to a leading order self-energy contribution to the 
quadratic part of the bosonic Hamiltonian 

S( <Z ) = -27Vr 2 ^ £fc+g /(£; fe ) 
fc 

Nr 2 

= — ^Wfo + aAfoJ ( 8 ) 

Fourier transforming to the momentum space, the 
bosonic part of the action reads 

ffbos = NY,A l {q)B^{q)A»{-q) 

qij,v 

without explicitly indicating the frequency dependence 
for the sake of simplicity and where (1,1/ = r, A, a. 




Fig 1. Leading-order self-energy contribution to the boson 
propagator from the four-leg vertex in H[ n t 



The matrix B /J, ' u , can be explicitely determined from 
Eqs.(||)-(||) and it is found that all elements are zero ex- 
cept for 
= rl [A° (1 - (l/2)/3 9 ) + aXl (1 - (l/2) 7 ,)] , B^ = 
B^ r = irlB a > a = (cjl + u>%)/cJa. 

The last ingredients of our perturbation theory are the 
vertices coupling the quasiparticles to the bosons 

A r (fc, q) = -2trl (e k+q/2 + £ fc _ g/2 ) (9) 

A x (k,q) = i (10) 
A a (k,q) = -2g (11) 

allowing to write the interaction part of the Hamiltonian 
in the following form 

H ™ = E C 1+|. AM ( fc > «) C fe-f <rW («) ■ ( 12 ) 
k.q.a 

The quasiparticle-boson interactions give rise to self- 
energy corrections to the boson propagators, which, at 
leading order in I/AT, are just fermionic bubbles with 
insertion of quasiparticle-boson vertices 

E k+ «-E k _«-iw m 
xA"(fc s9 ) A" (*,-«). (13) 

Once these self-energy corrections are taken into account, 
the boson propagator at leading order assumes the form 

D^(q,u, m ) = (A»(q,u> m )A u (-q,-u; m )) 

= N- 1 {2B + U{q,uj m ))-l (14) 

The factor 2 multiplying the boson matrix B arises from 
the fact that the bosonic fields in the presently used radial 
gauge are real. 

The above formal scheme allows to calculate the 
leading-order expressions of the effective scattering am- 
plitude both in the particle-hole channel 

T(k, k';q,w) = -J2^ (k', -«) (q, u) A" (k, q) . 

(15) 

and in the particle-particle channel 

x D» v {k- k',u)A v ^-*±*l,fc_ kfj (16) 

It should be noted that the boson propagators are of or- 
der 1 /N while the occurrence of a bare fermionic bubble 
leads to a spin summation and is therefore associated 
with a factor N. Thus, in this 1/N approach, the quasi- 
particle scattering amplitudes are residual interactions of 
order 1/N. 



4 



The form of the static density-density correlation func- 
tion at the leading order is 

a a' 

= P°( g ,« = 0) +5^X^(3^ = 0) 
x £>^(g,o; = 0)x° n fe^ = Q) 



where 



N ^ 

(7(7 



[n a {q)n a >{-q)) 



(17) 



(18) 



is the orbital bare density-density correlation function, 
and 

xJU«> w ) = ^ S<»a(g)X] C L A "( fe > 9)cfc+g«r')o (19) 



A simple inspection of the diagrammatic structure of the 
scattering amplitudes T and of the response functions 
P together with the observation that the static fermionic 
bubbles are non-singular function of q, allows to conclude 
that both quantities provide the same amount of infor- 
mation as far as the occurrence of the instabilities is con- 
cerned. In fact one can immediately recognize that, since 
a diverging response function can only arise from a di- 
verging boson propagator, also the scattering amplitude 
mediated by the same boson propagator would diverge at 
the same time. The same holds for the static Cooper scat- 
tering amplitude in the particle-particle channel, which, 
within our leading-order iM&f expansion coincides with 
the particle-hole amplitudeEl 



III. PHYSICAL PROPERTIES OF THE 
HUBBARD-HOLSTEIN MODEL 

A. Static properties 

The model introduced in the previous section is char- 
acterized by the contemporary presence of a very strong 
local interaction and a phonon-mediated attraction. The 
main point to be addressed here is the subtle interplay 
between the attractive and repulsive forces so as to clarify 
the origin of the effective interactions arising between the 
quasiparticles giving rise to instabilities for some values 
of the parameters. We will first analyze the static prop- 
erties, so that we first focus on the ui = limit of the 
effective interaction in the particle-hole channel Eq.(|l5|) 
between quasiparticles on the Fermi surface (fe = k Fl 
k = kp) . 

Since in general this quantity involves the calculation 
of the fermionic bubble s (|13| ) entering the expression of 
the boson propagator (|l4j) an explicit analytic evalua- 
tion of this scattering amplitude at finite momentum 



is not possible. The finite momentum analysis of the 
static scattering amplitude is then carried out numeri- 
cally. However we find instructive to present first the 
analytic results, which can be obtained in the small trans- 
ferred momentum limit (q — > 0) 



T q = lim lim T(kp, k' F ; q, to) 



q— >0 Ul— >o 




(20) 
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limit of the bubbles, IT^ = 



where the static q 

n M "(g — » 0,uj m = 0) can easily be evaluated by notic- 
ing that at T = 



W v (q^Q,u v 



°)=E 



df(E k ) 
dE k 



A M (fc, 0) A u (fc, 0) 



= -J2 5 (Ek~ E P ) A» (k, 0) A" (fe, 0) 

fc 

= -i/*A" {k F ,0)A u (k F ,0) 



(21) 



(22) 



where v* is the quasiparticle density of states at the 
Fermi level. In the last equality we use the property that 
the vertices are constant over the Fermi surface. Once 
the expressions for the bubbles are substituted in Eq. (|20| ) 
one obtains the final expression for the static scattering 
amplitude between the quasiparticles 



Nv*T n = 



Nv*Y u 



1 + Nv*T u 



Uv*e kF - A g 
1 + {Atv*e kF - A g ) 



(23) 



where we introduced the effective phonon-mediated e-e 
coupling A g = 2v*g 2 /uj$ In Eq.(p3|) the dynamical ef- 
fective scattering amplitude between the quasiparticles 
T u = T{kp, k,F] q = 0, a; — > 0) is identified by 



r. - i U„ - ± 



(24) 



This latter quantity represents the residual interaction 
between the quasiparticles on the Fermi surface when 
their mutual screening is not taken into .account. This 



screening effect is instead included in r g c3. The above 
expression for T u can also be easily obtained from the 
direct evaluation of Eq.(pfj|) by noticing that in the dy- 
namical limit the fermionic bubbles vanish identically. 
Thus, using = in Eq.(^i|) one again finds Nv*Y u] = 
4tv*e kF - Ag. 

At this point one recognizes that an instability can in 
principle take place when the Landau-Pomeranchuk cri- 
terion for the stability of a Fermi liquid Fq = Nv*Y UJ > 
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— 1 is violated leading to a negatively diverging total 
(static) scattering amplitude T q 



F° = Nv*T u = Atv*e kF - A- < -1 



(25) 



As a consequence a divergent compressibility is found 



On 



Nv* 



1 



> — 1 00. 



(26) 



which signals the occurrence of a PS. It is worth not- 
ing that the phonon parameters only enter the condition 
for PS via the combination \ g — 2v*g 2 /uq. This im- 
plies that, for a given electronic band structure, i.e. for a 
given i>*, the static instability is obtained for any phonon 
frequency provided g is suitably rescaled to keep the A g 
fixed (of course too large e-ph couplings put in jeopardy 
our weak coupling approach, where vertex correction be- 
yond the Migdal theorem are not included). 

It is worth remarking that the quantity £k F = 
(cos k X F + cosJbyf) + 2a (cos k X F cos k y p) is rather small 
at low dopingcll We carried out a detailed investigation 
of £k F as a function of the doping and of the hopping ratio 
a determining the most appropriate values of t' in order 
to reproduce within our single-band model Fermi surface 
shapes in reasonable agreement with those observed in 
some superconducting copper oxides. According to pre- 
vious analysescil, we found that the values of t' giving rise 
to reasonably shaped Fermi surfaces for 214 compounds 
are negative with t' ~ — 0.1672;. t' ~ —0.454 was instead 
adopted to reproduce Fermi surfaces in agreement with 
those observed in 123 and 2212 compounds. For these 
values of t 1 , £k F remains rather small in the low doping 
regime. Therefore, although a minumum critical value g c 
of the e-ph coupling is needed for the phase-separation 
instability to occur, its value can be small and it does 
not really provide a difficult condition to be satisfied in 
the real systemscl 

The main point to be stressed here is that, despite the 
infinite bare repulsion between the bare particles, a huge 
screening takes place in the system introducing attractive 
forces almost exactly balancing the repulsion and giving 
rise to a finite scattering amplitude between the quasi- 
particles. This residual effective interaction is repulsive 
in the absence of the e-ph coupling, whereas it may turn 
into an attraction when A g is large enough. In this situ- 
ation, which arises from the strongly interacting nature 
of the system, even this additional attraction can drive 
the system unstable. This finding matches well previous 
results obtained in different models and it emphasizes 
the general robustness of the PS concept in the context 
of strongly interacting systems. Moreover this finding 
supports the choice of the Hubbard-Holstein model as a 
simple paradigm to describe the physics of the PS in the 
context of strongly interacting systems. 

In principle the above analysis cannot exclude that an 
instability at a finite momentum can take place before 
PS, thus calling for a more general investigation. This 



latter has been carried out numerically through the ex- 
plicit calculation of both the static scattering amplitude 
at finite transferred momentum and the density-density 
response function. In Fig. 2 we report the full static 
scattering amplitude in the particle-particle channel. 
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Fig 2. (a) Doping dependence of the mass D in the static 
effective scattering amplitude close to the PS instability. The 
open squares in (a) indicate the values of D at various dop- 
ings and g phys = g cp h ys = 0.194/^ eV, t phys = 0.5eV, 
t' — — (l/6)t, and wophys = 0.04eV. For these parameters 
S c = 0.195. The solid line is a linear fit. (b) Static scatter- 
ing amplitude for the same parameters as in (a) as a function 
of the transferred momentum q in the (1,0) direction. The 
doping 5 = 0.2, 0.205, 0.22, 0.25, 0.28 increases from the lower 
solid line to the upper solid line. 

We choose this quantity because it will also enter the 
calculation of the Cooper instability reported below. We 
performed the calculation for various different values of 
the doping at a value of g larger than the minimum value 
required to have a PS instability. Here and throughout 
this paper we express the various quantities in physical 
units translated from the 1/N formalism with N = 2: 
tphys =t/N, t' phys =t'/N, g phys = g/VN, and w 0p h y s = 

UJQ. 

From the reported results it is natural to conclude that 
the divergent particle-particle scattering amplitude (i.e. 
a divergent boson propagator and a consequently diverg- 
ing density-density response function) occurs at a zero 
momentum transfer. According to Eq.(^6|) this diver- 
gency (r g — > — oo) leads to a diverging compressibility 
signaling a PS instability. 

As it can be seen in the g vs S phase diagram of Fig. 
3, for the t' = — 0.1674 case (a similar diagram is ob- 
tained when t' = — 0.454), we found that this behavior is 
rather generic at low and intermediate doping (continu- 
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d-wave channels. (Notice that A; > means attraction). 

The results are tabulated at various doping concen- 
trations for the case with t p h ys = 0.5eV, t — —1/61/, 
9phys = 0.192/V2eV, wophys = 0.04eV in Table I. 
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Extended s- and d-wave superconducting couplings for t p h ys = 
0.5eV, t = -l/6t', Sp hys = 0.192/ v / 2eV, u 0phy!s = 0.04eV at 
various dopings. The instability is at S c — 0.193. 



Fig 3. Phase diagram e-ph coupling g vs doping 5 t p h ys = 
0.5eV, t' = -(1/6)*, and w ph ys = 0.04eV. On the solid line 
the compressibility diverges, whereas the dotted line arises 
from the Maxwell construction. 

The instability line drawn in the phase diagram in- 
dicates where the static density-density response of the 
system becomes singular. A narrow dip in the instability 
solid line is due to the presence of a van Hove singular- 
ity enhancing z/*, thus favoring the instability. However, 
the system becomes unstable and phase separates before 
the solid line is reached, where k —> oo. A Maxwell con- 
struction (dashed line) is needed to determine the region 
where PS starts. Appendix A briefly describes the pro- 
cedure to carry out the Maxwell construction. 

We complete the static analysis of the model by in- 
vestigating the possibility of Cooper pairing. As already 
pointed out, in various models with strong e-e correla- 
tions superconductivity can appear in the proximity of a 
PS instability. This can easily be understood by looking 
at the large attraction arising in the particle-particle ef- 
fective scattering amplitude close to the instabilities (see 
Fig. 2). According to this simple observation and accord- 
ing to the previous experience in other strongly interact- 
ing models, we therefore investigated the Fermi surface 
average of the particle-particle scattering amplitude de- 
fined in Eq.@ 

Aj = - / dkS (E k - n) 9l {kf 

dkdk' [gi{k)T c (k,k';ui = 0)gi{k') 

xS(E k -ij.)S (E k , - (j,)] (27) 

with g Sl (k) — cos(k x ) + cos(k y ), and gd x {k) — cos(fc x ) — 
cos(fcy) projecting the interaction onto the s-wave and 
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With the set of parameters of Table I the critical dop- 
ing for the occurrence of the instability is 5 C = 0.192. 
Whereas the couplings Xd are found to be generally at- 
tractive near (and inside) the unstable region, s-wave 
Cooper instabilities are found only very close to the in- 
stability linet3. 

At first sight it may seem strange to investigate and 
to look for superconductivity in a phase-diagram region, 
which is made unaccessible by the Maxwell construction. 
However, one should remember that the phase diagram 
of Fig. 3 may be modified by several effects. First of 
all long-range Coulombic forces will spoil PS. This issue 
is the main topic of the present paper and will be dis- 
cussed in Section IV. Secondly temperature effects could 
restore the uniformity of the system. This issue cannot 
be reliably addressed within a slave-boson formalism and 
is beyond the scope of our work. Most importantly su- 
perconductivity itself can give rise to a more stable phase 
competing with and stabilizing PS. In this latter scenario 
a complex interplay between the PS instability and su- 
perconductivity will likely arise: The incipient instability 
originates superconductivity, which, in turn, prevents the 
instability to occurE-1. To substantiate these ideas a con- 
sistent analysis would be required of the feedback effects, 
which, however, only appear at higher order in the 1/N 
expansion. Work in this direction is in progress. 



A'; 



A e = 
q 1 + Nv*Tf, 



(30) 



where is the first term in the r.h.s. of Eq.(|24|) and 
AjJ, = 1. Eq.([30|) shows the difference between the dy- 
namic and static limits, which can be substantial for large 
v*. This is a generic feature of strongly correlated sys- 
tems, which one has to take into account since it can 
strongly affect the relevance of the e-ph coupling. 

The above small-g and small-w analysis can easily be 
extended to finite momenta and frequencies by an explicit 
numerical evaluation of the diagrams of Fig. 4. 




B. Dynamical properties 

The analysis carried out in the previous sections was 
purely static. From this analysis it turned out that the 
subtle balancing between repusive forces and attractive 
screening processes may have relevant physical effects. 
In this regard it immediately appears of great interest a 
dynamical investigation in order to fully clarify the way 
the various screening effects contribute to the physics of 
the system. _ 

From the standard Fermi-liquid theoryc3, one obtains 
the two relations connecting the density vertex A e (<7, u>) 
and the wavefunction renormalization z e in the dynamic 
and static limits 

z e A e (g = 0,w->0) = 1 (28) 

2 e A e (^0,^ = 0) = 1 - W r. (29) 

1 + F° {e) 

where = 2v*Tl J and r* are the Landau parameter 

and the dynamic (q — 0, u> — + 0) effective e-e scattering 
amplitude between the quasiparticles respectively, when 
only the e-e repulsion is taken into account. To explicitly 
keep memory of this limitation we append a suffix "e" to 
and to any quantity not involving phononic processes. 
In our specific model and within our leading-order 1 /N 
approximation one can recognize that the static and dy- 
namical e-ph vertices, A q and A w , are related by 



Fig 4. Leading-order in 1/N diagrammatic structure of 
the effective e-ph vertex dressed by electronic processes only: 
the wavy line is the slave-boson propagator only involving r 
and A bosons, the solid dot is the dressed e-ph vertex, the 
open dot is the bare e-ph vertex and the grey dots are the 
quasiparticle-slave-boson vertices. 

In these diagrams one considers the effects of the Sr.i 
and SXi boson fluctuations, keeping track at leading order 
in 1/N of the original infinite Hubbard repulsion U. More 
explicitly one has to evaluate 

A e -P h (fc,g;o; m ) = 

A a (k,q)-N ]T A»(k,q)D^(q,w m )lF' a (q,w m ) (31) 

ii,^— r, X 

For u> m — our slave-boson result is in perfect quanti- 
tativ e agr eement with the static analysis carried out in 
Refs.E3£J for the single-band infinite-U Hubbard model 
treated within a large- N expansion by means of Hubbard 
projectors. 

The results of a dynamical analysis of the e-ph ver- 
tex vs transferred Matsubara frequencies are reported in 
Figs. 5(a) and 5(b) for t' = -i/6 for a small [q = (0.2, 0)], 
and a sizable [q = (2.0, 0)] value of the transferred mo- 
mentum respectively. Similar results are obtained for 
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the t' = — 0.45i case. It is worth noting in Fig. 5(a) 
how rapidly the effective e-ph vertex increases as soon 
as the transferred frequency becomes larger than some 
screening scale cj scr of the order of v F q, with v* F oc 5. A 
much larger scale of the order of the bare bandwidth t is 
involved in the slower increase of the vertex at large mo- 
menta [Fig. 5(b)]. A detailed discussion of the screening 
scales is deferred to the next section, where this analysis 
will also be carried out in the presence of LRC forces. 




Fig 5. Effective e-ph vertex as a function of the transferred 
Matsubara frequency for t p h ys = 0.5eV, t' = — (l/6)t, and 
i^Ophys = 0.04eV. In (a) the transferred momentum is small, 
q = (0.2, 0), and it is large, q = (2.0, 0), in (b) 
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Fig 6. Effective scattering amplitude as a function of the 
transferred Matsubara frequency for t p h ys = 0.5eV, t' = 
-(l/6)t, ffphys = 0.205/\/2eV, and w 0p hy 8 = 0.04eV. Solid 
line: Small transferred momentum q — (0.15,0); dashed line: 
Large transferred momentum q — (1.75, 0). In the insets the 
attractive part r at tr = T — r rop (see text) are reported for 
both small (solid line) and large momenta (dashed line). 

Since dynamical effects strongly modify the behavior 
of the e-ph interaction, it seems natural to extend the dy- 
namical analysis to other relevant quatities of the system. 
In particular we investigated the finite frequency behav- 
ior of the effective Cooper scattering amplitudes between 
the quasiparticles on the Fermi surface r c (/s.p,Aj 
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The results are reported in Fig. 6 for the particle-particle 
scattering amplitudes as a function of the transferred 
Matsubara frequencies both for small and large momenta. 
In this case the bare e-ph coupling g p hy S = 0.194/v^eV 
and the doping 5 = 0.205 are tuned in order to place the 
system in the proximity of a q = instability occurring 

for ^cphys — 

0.194/^2 and S c = 0.195. For clarity we 
also represent in the inset the difference between the to- 
tal and the purely repulsive (i.e. involving the r and A 
bosons only) part of the scattering amplitude. For any 
momentum, this attractive part of the interaction lives 
on a frequency range of the order of the phononic energy. 
However, as expected, the large and the small momen- 
tum behavior are strikingly different. At small momenta 
the attraction mediated by the phonons is stronger be- 
cause the e-ph coupling is larger (cf. Fig. 5) and, close 
to the instability, it gives rise to a large attraction at low 
frequencies. This attraction is the vestige of the huge 
static attraction found close to the PS instability also in 
the particle-particle channel (cf. Fig. 2). On the other 
hand the attractive contribution in the large-momentum 
case (dotted curve in the inset of Fig. 6) is quit small 
for any frequency, because the effective coupling between 
the quasiparticles and the phonons is greatly reduced by 
the large-momentum screening [cf. Fig. 5(b)]. 

We finally like to comment on the behavior of the scat- 
tering amplitude for frequencies that are larger than both 
the phononic and electronic energy scales. In this case 
it can be easily checked that the fermionic polarization 
bubbles vanish as lu~ 2 thereby leading to a scattering am- 
plitude determined by the bare slave-boson propagators 
[the bare phonon propagator (B a ' a )~ 1 also vanishes for 
large frequencies] . A simple calculation based on Eq. ( |l6| ) 
and on Eq. (|l4|) with IL(q,u>) — ► gives the large-w sat- 
uration value of the scattering amplitude in the Cooper 
channel 



T (kF,k' F ;uj — > 00) 



1 

N 



■ite kf 



2rl 



(32) 



where Aq' 1 are defined in Eq.(|6|) and f3 q and j q are given 
after Eq.(||). A direct inspection of Fig. 6 shows that 
this saturation value is reached more rapidly in the small 
momentum transfer case, whereas a slower rise of T (w) 



9 



is found in the large momentum transfer case of Fig. 6. 
This result is a natural consequence of the disappear- 
ance of the phonons from the high frequency processes 
at (u> n > luq). Then only electronic processes determine 
the scattering like in the screening of the e-ph vertex de- 
scribed above [cf. Figs 5(a), (b)]. Again two energy scales 
appear to be relevant: A small u> scr ~ Sq at small mo- 
menta and a large uj scr ~ t at large momenta. These two 
scales set the frequency region above which the scattering 
tends to saturate. 

As far as the dynamic behavior of the scattering ampli- 
tude in the proximity of a PS instability, we also like, for 
completeness sake, to remind that a recent workQ also 
showed that, close to a q = instability, the effective 
dynamical scattering amplitude for real frequencies has 
a strongly singular behavior. In particular the effective 
interaction assumes the same form as the one obtained 
within the gauge-theory treatment of the t-J modeO 



1 



Bq 2 - iuj^ + D 



(33) 



with a mass D oc (5 — 6 C ), which vanishes linearly when, 
by varying the doping one approaches the q — in- 
stability line. U represents the almost momentum in- 
dependent repulsive contribution to T mediated by the 
r and A bosons. According to the spirit of the Fermi 
liquid theory, it may be interpreted as the residual re- 
pulsion surviving between the quasiparticles and arising 
from the infinite repulsion U between the bare electrons. 
Eq. (|33"|)establishes a connection between the presence of 
a q = charge instability and singular scattering, which 
could determine the anomalous normal state properties 
of the copper oxides. 

The analysis summarized in Figs. 5 and 6 is of ob- 
vious pertinence in a complete Eliashberg treatment of 
the superconductivity problem. In particular it is evi- 
dent that the huge enhancement of the attractive part 
of the scattering amplitude near the instability line can 
be responsible for large critical temperatures despite the 
small e-ph coupling. 

The closeness to a phase-separation instability appears 
therefore as a favorable condition in order to obtain high 
temperature superconductivity from a phonon-mediated 
attraction similarly to what suggested, in, the context of 
purely electronic pairing mechanismaL3El 



IV. THE EFFECT OF LONG-RANGE 
INTERACTIONS 

A. Formal extension 

Although we expect the effect of LRC forces to be 
most effective in the small transferred momentum case, 
where the underlying lattice structure is less visible, we 
cxplicitely kept into account the real symmetry of the 
square-lattice system. Therefore, to derive an explicit 



expression for the Coulombic potential in the spirit of 
the point-charge approximation we started from the dis- 
cretized form of the Laplace equation. Moreover, since we 
assume that our twodimcnsional model represents planes 
of a truely threedimensional lattice we also include a 
third spacial dimension. For clarity in this section we 
restore the explicit dependence of the lattice spacing a, 
which in the previous sections was set to unity in the 
square twodimensional lattice. In the third space direc- 
tion, instead, we assume the unit cell to have a lattice 
spacing d (i.e. we assume a tetragonal threedimensional 
lattice). In this scheme, the Laplace equation reads 



-ii £ 

ri=x,y 



J i—j-\-r) 



(34) 



d 2 



-e5(i-j) (35) 



where i,j are lattice sites and e± and eii are the dynamic 
dielectric constants perpendicularly and along the planes 
respectively. Fourier transforming one can easily obtain 



[cos(aq x ) + cos(aq y ) — 2] + cos(dq z ) — 1 



{a/df 



x^r0 q = -e (36) 

with e = e||/ej_, from which one gets the expression of the 
LRC potential in the threedimensional momentum space 

ed 2 _ x 
q = - o — [Mflx, %) + cos(q z d)} (37) 



where we defined 



Mlx,q v ) 



{a/df 



[cos{aq x ) 4- cos(aq y ) — 2] — 1 . (38) 



Since we are interested in the effects of the Coulomb po- 
tential on the square lattice planar system, we now trans- 
form from q z to the real space for the plane at z = 
obtaining 



6 q|| {z = 0) = - 



1 



2c 



± sj A 2 (q x ,q y ) - 1 



(39) 



Notice that this is the potential between electrons in a 
twodimensional lattice embedded in a threedimensional 
space and it diverges as q -1 for small transferred mo- 
menta, rather than q~ 2 as it happens for threedimen- 
sional electronic systems. This potential can be used in 
the Coulombic part of the Hamiltonian 



1 



where p q = J2k 
constant Vc 



^A 2 {q) - 1 



PqP-q 



(40) 



k+q a c k,<y an d the Coulombic coupling 
= e 2 d/(2ej_a 2 ). It should be noted that, 
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as it is customarily done, the sum does not include the 
zero-momentum component, since we are supposing that 
the diverging q — interaction between the electrons 
is canceled by the contribution of a uniform positively 
charged ionic background. Having in mind the supercon- 
ducting copper oxides of the 214 type, where d ~ 3a, 
iphys ~ 0.5eV, e|| w 30 and e± 5 one sees that Vc 
has to range from roughly 0.5eV to 3eV in order to have 
holes in neighboring CuC-2 cells repelling each other with 
a strength of 0.1-0.6eV. 

The Hamiltonian ( |40| ) can then be added to the Hamil- 
tonian of the Hubbard-Holstein model (^) and the prod- 
uct of four fermionic fields can be decoupled by means 
of a standard Hubbard-Stratonovich transformation. In 
this way one introduces a new real bosonic field Yi to be 
integrated over in the functional integral. Although this 
spatially fluctuating field does not have its own dynam- 
ics it acquires a frequency dependence via its coupling 
to the fermionic degrees of freedom. In particular, ex- 
tending the bosonic space A^(q) — (Sr q , SX q , a qi Y q ) (we 
dropped for simplicity the Matsubara frequency depen- 
dence of the bosonic fields) one can extend the formalism 
of Section II to include the effects of the Coulomb poten- 
tial represented by the Y boson. A direct calculation 
shows that the quasiparticle-boson vertex is 



A Y (k,q)=z 



(41) 



as expected since the Coulomb potential couples to the 
local electronic density in the same way as the boson A 
does. The bare boson propagator becomes a 4 x 4 matrix 
with an additional non-zero element: 



Y<Y = y/A*(q)-l 
2V C 



(42) 



We stress again that the Y boson does not have any uni- 
form q = component, which was discarded from the 
beginning in the sum of Eq.(j40|) and therefore it does 
not affect the mean-field results. 



B. Static properties 

The introduction of a LRC potential greatly affects 
the coupling between phonons and electrons particu- 
larly for small momenta, where the Coulombic forces are 
most effective. In fact, since in the Hubbard-Holstein 
model phonons couple to the local electron densities 
like the Coulombic potential, the e-ph coupling is effec- 
tively screened by particle-hole pairs created by the e-e 
Coulombic scattering processes and the effective e-ph ver- 
tex is largely suppressed in the small- momentum limit. 
This effect appears in a direct calculation of the expres- 
sion (|l]) for the effective e-ph vertex in the presence of 



LRC forces, A^ fi ph (fc, q; ui m ). Now fi, if = r, X, Y in (|3lj) 



so as to include the screening due to the fluctuations of 
the "long-range" boson Y. In particular Fig. 7 shows 



the momentum dependence of the static (oj n = 0) e-ph 
vertex. The comparison with the same quantity in the 
absence of LRC forces shows a striking difference in the 
low-momentum region, where the Coulombic screening 
leads to a vanishing e-ph vertex. 
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Fig 7. Static effective e-ph vertex as a function of the trans- 
ferred momentum in the (1,0) direction for i p h ys = 0.5eV, 
t' = — (l/6)t, <x>o P hys = 0.04eV. The dotted line is in the ab- 
sence of LRC forces (Vbphys = 0); the solid line is in the 
presence of LRC forces with Vbphys = 0.55eV 

This difference can easily be understood by a three- 
step evaluation of A c £~ I ^ h (k,q;u! m ). One can first go 
through the calculation of the vertex in the absence of 
LRC forces [Eq.©]. Then one can introduce a density- 
density fermionic bubble dressed by all electronic short- 
range processes 

Ar 2 P°(q,to m ) 



Il eSR (q,u, r , 



(43) 



dot £ 2 >< 2 (q, CJ m ) 

where D 2x2 is the 2x2 sector of the boson propaga- 
tor only including the r and A bosons, representing the 
purely electronic short-range processes. Finally one per- 
forms a resummation only including Y-boson fluctuations 
dressed by the Yl eSR density-density fermionic bubbles. 
The result is 

A c -P h (fc, o;Wm ) 



A Li? ph (M; w ™) 



l+n e ^(o;w m )- 



(44) 



s/A*{q)-l 

This expression readily shows the suppression of the 
short-range-only e-ph vertex appearing in the numera- 
tor and clearly displays the vanishing of the LR ver- 
tex for q — * 0, when the Coulombic potential diverges 
(A2(o)-l)-i^|q|-i ; 

The important physical consequence of the above sup- 
pression of the effective e-ph vertex is that in the pres- 
ence of LRC forces the long- wavelength density fluctua- 
tions are decoupled from the phonons and become unable 
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to drive a low-momentum instability. Nevertheless the 
possibility of finite-momentum instabilities still remains 
open and a detailed reanalysis of the phase diagram in 
the presence of the LRC potential is needed. This anal- 
ysis can be carried out by direct numerical evaluation 
of the static density-density response function and of its 
divergencies. A direct inspection to the diagrammatic 
structure of the leading-order in 1 /N density-density re- 
sponse function shows that one can include the Coulom- 
bic effects both by extending the fi, v summations in Eq. 
( |l7| ) to include the Y boson or, alternatively, by first cal- 
culating the short-range response function P SR by only 
including the r, A and a bosons and then to resum with 
the bare (YY) propagator. The result is then 



P LR (q,LJ = 0) 



P bH (q,u = 0) 



v c 



=P SR (q,U = 0) 



(45) 



From this expression one can easily see that a positively 
diverging P SR (q,ui = 0) no longer gives a diverging LR 
densit y-density r esponse function. In particular, since 
Vc/(\/ A 2 (q) — 1) — * oo for q — > 0, the compressibility 
always vanishes as it should in a Coulomb gas and a 
PS instability is now ruled out. However, some finite-^ 
instabilities are still possible in the system when 



P SR (q,uj = 0) = -N(^/AZ(q)-l)/V c 



(46) 



leading to a divergent P LR . This is possible in principle, 
since inside the PS region for the model without LRC 
forces, P SR has a simple pole. Therefore P SR has a neg- 
ative branch from zero momentum up to sizable momenta 
(if the parameters are choosen to be inside enough into 
the unstable region) and the condition ( f46| ) can be satis- 
fied. At this point, then, an instability occurs despite the 
stabilizing effect of the LRC potential and an incommen- 
surate charge-density-wave (CDW) phase takes place in 
the system. 

The g vs 8 phase diagram for various values of Vc is 
reported in Fig. 8, where the continuous line represents 
the place at which P LR (q = q c ,uj — 0) diverges separat- 
ing the stable uniform region from the unstable region 
where an incommensurate CDW is expected to form. It 
is important to notice that, when the settling of a CDW 
phase is given by a second order quantum transition, no 
Maxwell construction is needed to determine the stability 
region. 
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Fig 8. Phase diagram e-ph coupling g vs doping 8 with 
iphys = 0.5eV, t! = — (l/6)i, wophys = 0.04eV, and in the 
presence of LRC forces with Vcphys = 0.55eV (lowest curve), 
Lbphys = 1.65eV (middle curve) and Vc P hy S = 3.3eV (upper 
curve) . 

As it can be seen the effect of LRC forces is stronger 
at low doping. This is so because in this region the poles 
of P SR tend to occur at low momenta and are more ef- 
fectively stabilized by the Coulombic potential. 

We found that the momenta q c at which the divergen- 
cies in P LR occur obviously depend on the point g vs 8 
and on the strength of the Coulomb force Vc, but are 
generically sizable indicating that the wavelength of the 
expected CDW phase is of a few unit cells. 

We performed an extensive analysis of the finite-g in- 
stabilities for various values of doping, f'/f, g, and Vc- 
It turns out that, for reasonable values of these parame- 
ters, the instabilities always occur at or close to the (1,0) 
and (0,1) directions. This effect is due to the momen- 
tum structure of the short-range density-density response 
function, which is enhanced by the large density of states 
in the (1,0) and (0,1) directions. 

It is interesting to notice that, close to the finite-g in- 
stability, the static effective scattering amplitude in the 
particle-hole as well as in the particle-particle channel 
diverge at finite momentum transfer likeo 



r(^ = o) = -— 



D> + B'\q 



(47) 



with D' oc (8 — 8 C ) being a mass linearly vanishing by 
approaching S c . The evolution of T (Uf ,k' F ,uj — 0) in 



approaching the unstable region by doping variations is 
displayed in Fig. 9. 
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111(a), a static Cooper instability takes place when the 
Fermi-surface average of — T c (kp, k' F , uj — 0) becomes 
positive [cf Eq.(p7|)] and this seems likely to occur close 
to finite-momentum instabilities. The resulting coupling 
constants are reported in Table II. 
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Extended s- and d-wave superconducting couplings for t p h ys = 
0.5eV, t = -1/61/, Pphy S = 0.260/\/2eV, cJo p hy a = 0.04eV and 
Vbphys = 0.55eV at various dopings. The instability is at 
8 C = 0.299. 



Fig 9. (a) Mass D' of the effective static scattering am- 
plitude as a function of S — S c for t p h ys = 0.5eV t' = — l/6i, 
Vcphys = 0.55eV, cu 0phys = 0.04eV and g phya = 0.240eV. The 
open squares indicate the values of D at various dopings and 
the streight line is a linear fit. (b) Static scattering ampli- 
tude for the same parameters as in (a) as a function of the 
transferred momentum q in the q c w (±0.28/a, ±0.86/a), di- 
rection. The doping S — 0.195,0.2,0.205,0.22 increases from 
the lower solid line to the upper dotted line. 

Fig. 10 displays T c (kF, k' F , uj = 0) over a large portion 
of the Brillouin zone for parameter values close to those of 
a finitc-q instability. It is evident that a deep attractive 
interaction between the quasiparticles arises on a broad 
momentum region. 




Fig 10. Momentum dependence of the static scattering 
amplitude for the same parameters as in Fig. 9 at 8 — 0.195. 

This fact is particularly remarkable in relation to the 
occurrence of Cooper instabilities. As seen in Subsection 
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As expected, positive (i.e. attractive) coupling con- 
stants arise close to the incommensurate CDW instabil- 
ity once more supporting the idea that high-temperature 
superconductivity could arise in the proximity of a charge 
instability. It is also worth emphasizing that, like in the 
proximity of PS, such a pairing instability takes place 
as an effect of a q-independent phonon-induced attrac- 
tion, which, close to the CDW instability becomes highly 
structured in momentum space, thereby opening the way 
to pairing symmetries other than the simple s-wave. This 
fact accounts for the attractive couplings also found in 
the d-wave channel. In particular it turns out that the 
rf-wave symmetry of the order parameter is able to take 
good advantage of the strong small-q attraction and of 
the local (large-g) repulsion. 

In the region of small-intermediate <?'s and dopings our 
analysis indicates the sure existence of d-wave pairing in 
sizable regions near the instability, whereas the occur- 
rence of s-wave pairing takes place in a much narrower 
region. However, it should be emphasized that the pres- 
ence of a s-wave static Cooper instability only in a nar- 
row region, by no means excludes the possibility of hav- 
ing s-wave superconductivity in a much larger area of 
our phase diagram. An appropriate Eliashberg dynami- 
cal analysis would be required to draw a firm conclusion, 
specially in the light of the results reported in the next 
section showing a strong frequency dependence of the ef- 
fective interaction between the quasiparticles. Of course 
the same applies to the attraction in the d-wave channels, 
which could also be greatly favored by dynamical effects. 

C. Dynamical properties 

We carried out a dynamical analysis of the effective 
e-ph vertex in the presence of LRC forces. The dynam- 
ical behavior of the e-ph vertex in the presence of LRC 
forces is not very different from the behavior observed 
for Vc = 0. In both cases two different screening regimes 
take place for small and large momenta. Specifically, 
when LRC forces are absent, we found that at small mo- 
menta the screening energy uj scr is proportional to the 
doping and depends linearly on the exchanged momen- 
tum uj scr 1=3 v* F q. Instead, in the presence of LRC forces, 
it is found that 

User ~ ( 48 ) 

The linear momentum dependence in the short-range- 
only case is compatible with screening processes associ- 
ated to both the particle-hole continuum and the zero 
sound. In the presence of LRC forces, the square-root 
momentum dependence is a clear indication that the plas- 
mon collective modecZl sets the cut-off energy at which 
small-momentum screening processes cease to be rele- 
vant. 

On the other hand we found that at large momenta, 
a different energy scale rules the screening processes, 



which, both with and without LRC forces, is of the order 
of the bare electron hopping t uj scr w jt. This indicates 
that the local (large-g) physics is governed by incoherent 
bare-electron processes with typical energies of order t 
much larger than the typical energies of order v* F q ~ tS 
ruling the coherent quasiparticle processes. 

The screening effects leading to different behaviors of 
the effective e-ph vertex at small and large momenta also 
affect the dynamical scattering amplitude. This generic 
effect is more evident in the proximity of an instability. 
Therefore in Fig. 11, we show the numerical evaluation of 
the dynamical scattering amplitude. A large attraction 
at low frequency is found for transferred momenta close 
to q c as a remnant of the infinite static attraction taking 
place at the instability. 
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Fig 11. Effective scattering amplitude as a function of 
the transferred Matsubara frequency for t p h ys = 0.5eV, 
t' = -(l/6)i, ffphy S = 0.240/v^eV, U) 0phys = 0.04eV, and 
Vcphys = 0.55eV. The doping 5 = 0.2 is close to the crit- 
ical value 5 C = 0.195. Solid line: Total effective scattering 
amplitude V; dashed line: Repulsive (slave- and Coulomb 
force bosons only) effective scattering amplitude r rop ; dot- 
ted lines: Attractive part of the scattering amplitude r a ttr = 
r — r rcp . The transferred momenta are in the instabil- 
ity direction (0.28,0.86). (a) Sizable transferred momentum 
\q\ = \q c \ — 0.9; (b) small transferred momentum \q\ = 0.15. 
The insets are enlargements of the attractive parts of the ef- 
fective scattering amplitudes. 

This attraction is, however, rapidly spoiled by increas- 
ing the Matsubara frequency. This is so because at fi- 
nite q « q c the dynamical screening of the e-ph ver- 
tex strongly suppresses the phononic attractive part of 
the scattering amplitude. Therefore, as soon as the fre- 
quency spoils the static attraction due to the instability, 
the attractive part of T c (dotted line in Figs. 11) rapidly 
vanishes. An opposite behavior occurs at small momenta, 
where no instability occurs and therefore no static attrac- 
tion is present. The usual e-ph attraction only appears 
at a finite frequency around loq > to scr ~ qS, where the 
e-ph vertex is not severely suppressed. 

We conclude this subsection by recalling the result of 
Rcf.El where a real-frequency analysis was carried out 
finding strong singular scattering between the quasipar- 
ticles of the form 

T(q,io)^U — (49) 

Wq — VjJ 

where w q = D' + B'|q - q c | 2 with D' oU8 - 5 C ). This 
scattering is of the form proposed in Ref.El to explain the 
anomalous normal state properties of the superconduct- 
ing copper oxides. Some substantial differences are, how- 
ever, worth being emphasized. First of all the origin of 
the singular behavior is not necessarily related to a mag- 
netic scattering mechanism, in so far it arises from the 
closeness to a charge instability. The mechanism driving 
the instability can be of various nature, magnetic, exci- 
tonic or, as in the present model, phononic. Secondly, as 
it can be seen in Fig. 10, a much more isotropic region 
of large scattering arises in the present context thereby 
bypassing the objection rised in RefO that only a few 
"hot" points on the Fermi surface undergo such strong 
scattering processes. 

D. The collective modes 

The analysis of the dynamical behavior of the model 
can be completed by investigating the collective modes 
present in the system. This study is particularly relevant 
for a deep understanding of the dynamical mechanisms 



ruling the instability formation. To emphasize similari- 
ties and differences we consider both the case with and 
without the LRC potential. 

It was repeatedly pointed out in the context of mod- 
els with short-range interactions only, that PS occurs 
without a softening of a massive collective mode. This 
was first established for the charge-transfer mode in a 
three-band Hubbard model with NN Xjoulombic repul- 
sion between copper and oxygen holesliS and it was also 
confirmed in the case of a three-band Hubbard-Holstein 
model with holes coupled to an optical phononEJ. 
A real-frequency analysis of both the imaginary part 
of the dynamical density-density correlation function 
Im(nn)(q,u) and of the poles of the phonon propagator 
D a,a (q,ui) shows that also in the single-band infinite-U 
Hubbard-Holstein model, PS occurs due to the phonon- 
mediated attraction between quasiparticles, which pulls 
the zero-sound mode into the particle-hole continuum. 
When the attraction is large enough, the velocity of this 
strongly damped mode vanishes eventually driving the 
system unstable at zero momentum. On the contrary, the 
phonon frequency, although sizably renormalized by the 
e-ph interaction stays finite all over the Brillouin zone. 

The scenario is strongly modified in the presence of 
LRC forces, when the instability takes place at finite 
momenta. For such finite momenta, in the absence of 
e-ph interaction the phonon is the mode at lower energy. 
With the introduction of a finite e-ph coupling, it turns 
out that the phonon is somewhat softened. In partic- 
ular the softening is complete, i.e. the phonon energy 
vanishes when the CDW instability takes place. Again a 
comparative real-frequency analysis of Im(nn)(q,ui) and 
of ImD a ' a (q, ui) allows to identify the nature of the mode 
entering the continuum and producing the large enhance- 
ment of the absorption at low-frequency. Fig. 12(a) re- 
ports the behavior of the phonon frequency [extracted 
from ImD 44 (q,u>)] as a function of momenta in the spe- 
cific direction at which the CDW instability takes place 
[(±0.2, ±0.8) and (±0.8, ±0.2) for the parameters related 
to 214 systems] at and above the critical doping value. 
In Fig. 12(b) the phonon dispersion are reported for the 
(±1,±1) direction for the same values of doping. 
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Fig 12. Phonon dispersion curves (a) in the instability 
(0.28,0.86) direction and (b) in the (1,1) direction, for t p h ys = 
0.5eV, t' = -(l/6)t, <7 P hy S = 0.240/V/2CV, w ph ys = 0.04eV, 
and Vcphys = 0.55eV. The solid curves correspond to the crit- 
ical doping 8 C = 0.195, the dot-dashed and the dotted curves 
correspond to 8 — 0.21 and 8 = 0.3 respectively. 

Far from the critical momentum the phonon frequency 
is sizably reduced with respect to its bare value. How- 
ever, only a minor doping dependence of this substantial 
reduction is visible in generic positions in the Brillouin 
zone making the observation of the softening effect rather 
difficult. On the other hand, a strong doping dependence 
of the phonon dispersion is apparent in the proximity of 
the critical momentum. In this momentum region the 
phonon completely softens at 8 = 6 C = 0.195. By increas- 
ing the doping, i.e. by moving away from the instability 
the softening becomes incomplete and progressively less 
important although it seems to stay substantial up to 
S = 0.3. 

Some considerations are in order on the above anal- 
ysis of the collective modes in the proximity of charge 
instabilitites. First of all a clear distinction can be made 
between the behavior of the phonon mode in the absence 
and in the presence of LRC forces. At a PS instabil- 
ity, the phonon always stays massive and only mixes 
strongly with a low energy zero sound mode at finite 
but very small momenta (q < u>o/vp). On the contrary 
the phonon becomes completely soft at the CDW insta- 
bility. Therefore it would be natural to indicate neutron- 
scattering experiments in copper oxides as crucial tests 
in order to determine how close these material are to a 
phonon-mediated CDW instability: if a substantial soft- 
ening of some phonon mode is found at some momenta, 
this would be a clear indication of a frustrated PS almost 
leading to CDW. Unfortunately, our schematic single- 
band Hubbard-Holstein model is not detailed enough to 



provide indications of the modes, which could undergo a 
visible softening in real materials. Moreover, the strong 
reduction of the phonon frequency is substantial on a siz- 
able momentum range, but is only found for doping close 
to the critical 5 C , and this fact could render the search for 
such an effect quite a difficult task. An additional diffi- 
culty is that, due to the strong mixing of electronic and 
phononic degrees of freedom, the softened phonon modes 
would be greatly broadened and they could well be seen 
only as an increase of spectral weight at low frequen- 
cies, quite differently from the theoretical 1/N picture 
reported in Fig. 12. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we analyzed the screening processes and 
the occurrence of instabilities in the Hubbard-Holstein 
model in the framework of Fermi-liquid theory. In partic- 
ular we investigated the role of LRC forces in stabilizing 
PS and producing incommensurate CDW instabilities. 

However, some limitations should be kept clear for a 
correct understanding of the scenario here presented. We 
carried out a leading order analysis in 1/N, which defi- 
nitely neglects some effects, which are relevant in a com- 
plete quantitative understanding of the real materials. 
First of all our approach is designed to deal with charge 
degrees of freedom, but lacks a correct treatment of the 
spin degrees of freedom, which turn out to be crucial in 
the low-doping phase of the copper oxides. Therefore 
antiferromagnetic correlations and the interplay between 
spin and phonon degrees of freedom are absent in our 
large-N, slave-boson model. These effects only appear at 
higher order in 1/N. Secondly, our leading-order expan- 
sion does not allow for phonon vertex corrections beyond 
the Migdal theorem, although in the low-doping region 
the quasiparticlc bands become very narrow thus leading 
to a violation of the condition Ep 3> u>o- The absence 
of these vertex corrections obviously rules out the possi- 
bility of a correct observation of multi-phonon polaronic 
effects. 

In the absence of LRC forces, the scenario found here 
is consistent with previous results obtained in the three- 
band Hubbard-Holstein model, where PS was also found. 
The main features observed here are i) a generic suppres- 
sion of the e-ph vertex due to electronic screening also 
resulting in a strong dependence on the vpq/co ratio; ii) 
the persistent possibility of PS arising from the phonon- 
mediated attraction; iii) a singular behaviour of the effec- 
tive interaction between quasiparticles at low momenta 
close to the PS region. 

The main achievement of the present work is, however, 
the analysis of the model in the presence of LRC forces. 
In particular it is remarkable that, when LRC forces are 
included in the model, PS is spoiled, but finite momen- 
tum instabilities still take place on substantial and quite 
physical regions of the parameter space. Singular effec- 
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tive interactions are again obtained at finite momenta in 
the proximity of the CDW instability; iv) Cooper pairing 
both in the s- and the (i-wave channels is present already 
in the static limit for systems close enough to the unsta- 
ble regions, both in the presence and in the absence of 
LRC forces. 

The present model provides a rather simple playground 
to investigate the generic properties of electronic systems 
close to charge instabilities. Indeed, while some prop- 
erties like, e.g., the behavior of the phonon mode, are 
strictly related to the phononic nature of the interac- 
tion here considered, the main results obtained here are 
generic of modelsEZl showing charge instabilities. In par- 
ticular we believe that the strict relation between charge 
instabilities, strong scattering and Cooper pair formation 
[points iii) and iv) above] is a generic feature of strongly 
correlated electron systems irrespective of the underlying 
physical mechanisms leading to PS or to incommensurate 
CDW. 
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APPENDIX A: THE MAXWELL 
CONSTRUCTION 

In order to perform the Maxwell construction, we mod- 
ify the model in Eq. ([!]) coupling the phonons to the 
full electron density rather than to the density fluctua- 
tions. In this way an e-ph coupling is effective already 
at mean-field level, where a non-zero mean-field value of 



in arises. The mean-field 



the phonon field ag = (a) 
Hamiltonian acquires the form 



H 'mf = Y] E kc\ a c ka - (no - A + 2ga Q ) V c\ a c k 



+ iVA ( - - ) + NwoaZ 



(Al) 



Minimizing the mean-field free energy with respect to ao 
allows to obtain the self-consistency equation for ao 



Nco ao = NgJ2f (E(k)) = g- (1 - S) 



(A2) 



This determines the phonon- induced shift of the chemical 
potential giving rise to a doping dependent correction to 
this quantity. Owing to the standard expression for the 
compressibility K = I s = — the q = instability can 
be read directly from the /i vs n curves: a stationary 
point in /j,(5) corresponds to an infinite compressibility 
and a standard Maxwell construction in [i(5) determines 
the region where the system is in a single phase. 
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